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MODELING  SWITCHED  RELUCTANCE  MOTORS 
UNDER  MULTI-PHASE  EXCITATION 


INTRODUCTION 


The  Second  Annual  Report  has  been  divided  into  three  parts  in  order  to  summarize  the 
activities  performed  on  this  ONR  grant  from  May  1,  1999,  until  May  31,  2000.  The  three  parts 
are  the  following: 

•  Part  1  entitled  “Modeling  Switched  Reluctance  Motors  under  Multi-Phase  Excitation”  covers 
“modeling”  activities  performed  under  Task  (b)  -  “Design  of  the  SRM”  -  of  this  ONR  Grant. 

•  Part  2  entitled  “An  Enhanced  Simple  Method  for  Designing  Switched  Reluctance  Motors 
under  Multi-Phase  Excitation”  covers  “design”  activities  performed  under  Task  (b)  -  “Design 
of  the  SRM”  -  of  this  ONR  Grant. 

•  Part  3  entitled  “Modeling  Switched  Reluctance  Motors  under  Multi-Phase  Excitation”  covers 
“control”  activities  performed  under  Task  (c)  -  “Design  of  the  SRM  Converter”  -  of  this 
ONR  Grant. 

The  reminder  of  this  report  refers  to  Part  1  of  the  Second  Annual  Report.  This  Part  1  has 
been  divided  into  the  following  sections: 

•  Section  A  evaluates  different  modeling  techniques  available  for  Switched  Reluctance  Motors 
(SRM)  operating  under  single-phase  excitation. 

•  Section  B  analyzes  two  specific  techniques  for  single-phase  excitation;  namely,  the  gage 
curve  and  ‘look-up-table’  modeling  techniques. 

•  Section  C  addresses  two  specific  techniques  for  multi-phase  excitation  of  SRM;  one 
technique  is  based  on  the  ‘look-up-table’  used  for  single-phase  excitation  and  the  other 
implements  an  Artificial  Neural  Network  (ANN)  to  predict  the  flux  linkage  and 
electromagnetic  torque. 


A  MODELING  TECHNIQUES  FOR  SINGLE-PHASE- 
EXCITED  SRM 


Simulation  must  be  incorporated  as  an  integral  part  of  the  design  procedure  of  a  SRM 
since  there  are  no  general  models  or  equations  describing  the  SRM  dynamic  behavior.  Moreover, 
SRM  does  not  have  the  'steady-state'  operation  found  in  traditional  electric  machines.  The  state 
variables  (i.e.  phase  currents  and  fluxes)  never  reach  their  constant  values  since  its  principle  of 
operation  is  based  on  a  series  of  'strokes'. 

In  the  design  of  traditional  electric  machines,  static  and  dynamic  simulations  have  been 
executed  as  separated  tasks  usually  performed  by  different  engineers.  In  the  case  of  the  SRM, 
however,  both  types  of  simulation  should  be  combined  together  since  the  SRM  drive 
(motor/converter  combination)  has  to  be  optimized  for  a  particular  application  (torque  vs.  speed 
characteristic)  to  achieve  high-performance  operation  [A.  1]. 

The  SRM  optimization  procedure  is  an  iterative  process  among  mathematical 
expressions,  static  and  dynamic  simulations  as  shown  in  Figure  A-l.  Therefore,  it  is  desirable  to 
use  a  Computer-Aided  Design  (CAD)  tool  to  evaluate  the  behavior  of  different  design  parameters 
obtained  from  design  equations  and  FEA  simulation  as  well  as  different  converter  topologies  and 
control  strategies.  This  CAD  tool  has  to  predict  with  high  accuracy  the  relevant  performance 
parameters  such  as  the  time  waveforms  of  flux  linkage  X(t),  phase  current  i(t),  rotor  speed  co(t), 
rotor  position  Q(t)  and  torque  T(t),  from  a  reduced  and  representative  group  of  known  data  or 
measurements.  Furthermore,  the  selected  CAD  tool  should  be  computationally  efficient  to  reduce 
the  time  required  for  the  optimization  process. 

This  Section  introduces  the  basic  nomenclature  and  mathematical  description  of  the 
single-phase  excited  SRM.  Next,  it  presents  and  classifies  different  dynamic  simulation 
techniques.  Finally,  it  suggests  the  SRM  modeling  technique  and  simulation  package  based  on  the 
simulation  goals. 


SRM  Specifications 


Figure  A-l  SRM  optimization  process. 


Figure  A-2  shows  the  basic  electromagnetic  structure  for  a  6/4  SRM.  Each  stator  pole  has 
a  coil  wound;  coils  are  connected  in  such  way  that  their  magnetic  fluxes  are  added  to  form  a 
phase  of  the  SRM.  Usually,  coils  of  opposite  poles  are  connected  in  series  to  form  a  phase  as 
shown  in  Figure  A-2.  A  phase  is  at  the  aligned  position  when  any  pair  of  rotor  poles  is  exactly 
aligned  with  respect  to  the  excited  stator  poles  (Phase  1  in  Figure  A-2).  At  alignment,  the 
magnetic  reluctance  is  at  its  minimum  saturating  the  flux  path  at  lower  current  levels;  therefore, 
the  inductance  is  at  its  maximum.  A  phase  is  at  the  maligned  position  when  the  interpolar  axis  of 
any  pair  of  rotor  poles  is  exactly  aligned  with  a  line  located  at  midway  between  two  adjacent 
stator  poles  (Phase  2  in  Figure  A-2).  It  is  noted  that  the  inductance  value  of  a  phase  varies  widely 
as  a  function  of  both  rotor  position  and  phase  current. 


Figure  A-2  Basic  SRM  electromagnetic  structure. 


The  SRM  works  in  a  sequence  of  'strokes'  where  the  angle  is  given  by  the  following 
equation: 


6 


ST  ~ 


360° 

mNR 


(A-l) 


where  m  is  the  number  of  phases  and  Nr  is  the  number  of  rotor  poles. 

Typically,  the  number  of  stator  poles  exceeds  the  number  of  rotor  poles.  The  torque- 
production  capability  of  one  rotor  pole  is  given  by 


The  torque  capability  angle  is  divided  into  two  regions.  In  the  first  region,  the  SRM  has 
positive  torque  production  (pM  motoring  torque).  In  the  second  region,  the  SRM  has  negative 
torque  production  (pG  generating  torque).  This  is  illustrated  on  Figure  A-3.  In  the  case  of  a  6/4 
SRM,  P  =  9(f.  Therefore,  in  the  first  45°  (-45°  <  Pm  <  0°)  the  SRM  develops  a  positive  torque,  and 
in  the  second  45°  (0°  <  pG  ^  45°)  the  SRM  develops  a  negative  torque. 


Figure  A-3  Stroke  angle  and  torque  capability  angle  for  a  6/4  SRM. 


The  SRM  works  under  single-phase  excitation  if  each  stator  pole  is  excited  within  the 
stroke  angle.  However,  normally,  Pm  =  >dsr  given  room  for  multi -phase  excitation.  If  the  stator 
pole  is  excited  an  angle  larger  than  the  stroke  angle,  the  SRM  works  under  multi-phase  excitation. 
In  the  case  of  the  6/4  SRM  presented  in  Figure  A-3,  dsr  =  30°,  P  =  90°  and  pM  =  45°;  therefore, 
two  phases  can  be  excited  during  15°  and  thus  developing  positive  torque. 

If  current  is  applied  to  a  phase  when  the  rotor  is  not  at  the  aligned  position,  the  rotor  will 
be  attracted  to  the  aligned  position  developing  a  torque  ideally  proportional  to  the  square  of  the 
phase  current  and  the  inductance  rate  of  change  with  respect  to  rotor  position.  The  most  general 
expression  for  the  torque  produced  by  one  phase  at  any  rotor  position  is  given  by: 

(A-3) 

i=const 

where  W'  is  the  coenergy ,  i  is  the  phase  current,  and  6  is  the  rotor  position.  At  any  position,  the 
coenergy  is  the  area  below  the  magnetization  curve  Mi,  0)  as  shown  in  Figure  A-4.  Consequently, 
the  coenergy  is  expressed  as: 
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where  X  is  the  flux-linkage. 
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A-4  Instantaneous  torque  as  the  rate  of  change  of  coenergy  at  a  constant  current. 


During  the  rotor  displacement  Ad  =62  -  6,,  there  is  an  exchange  of  energy,  AW,  between 
the  supplied  energy  and  the  stored  field  energy.  To  insure  that  the  mechanical  work  done  is  equal 
to  the  change  in  coenergy,  constant  current  is  assumed  as  the  rotor  moves  through  Ad.  The  torque 
can  be  thought  as  the  work  AW  divided  by  rotor  displacement  Ad.  Finally,  it  must  be  noted  that 
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where  co  is  the  SRM  rotor  speed. 


Most  SRM,  even  those  ones  operating  under  single-phase  excitation  such  as  the  6/4 
SRM,  operate  under  multi-phase  conditions  for  small  periods  of  time.  An  accurate  SRM 
modeling  technique  should  take  into  account  multi-phase  operation.  However,  single-phase 
excitation  modeling  techniques  are  still  being  used  due  to  several  reasons: 

•  Single-phase  modeling  techniques  are  simpler  and  faster  than  multi-phase  excitation 
modeling  techniques. 

•  The  number  of  data  points  required  for  multi-phase  modeling  technique  increases 
exponentially  with  the  number  of  excited  phases. 

•  Multi-phase  operation  can  be  modeled  as  a  superposition  of  several  single-phases 
assuming  linearity. 


It  is  not  the  case  of  the  SRM  since  it  works  under  high  magnetic  saturation;  therefore,  its 
description  is  non-linear.  However,  this  type  of  approximation  frequently  provides  enough 
accuracy  as  a  first  approach  to  model  the  SRM  dynamic  behavior  or  the  control  technique.  For 
higher  modeling  accuracy,  specific  multi-phase  excitation  modeling  techniques  are  required. 

Most  of  the  simulation  techniques  use  a  curve-fitting  technique  to  represent  the  flux- 
linkage  surface.  Such  a  technique  requires  a  large  number  of  measured  or  calculated  data  points. 
Figure  A-5  shows  the  flux-linkage  characteristic  of  a  6/4  SRM  calculated  using  FEA  simulations. 
The  FEA  of  a  SRM  is  a  complex  and  time-consuming  but  unavoidable  task  in  the  design  of  high- 
performance  SRM  drives.  The  code  writing  and  debugging  (setup  time)  of  a  FEA  simulation 
takes  approximately  a  week  (40  hours).  A  single  simulation  point  in  a  SUN'8'  workstation  machine 
(100  MHz  clock  CPU  and  120  Mbits  RAM)  takes  approximately  4  minutes. 

Figure  A-5  has  225  points  (15  different  rotor  positions  times  15  different  current  levels). 
Consequently,  once  that  the  FEA  code  is  debugged,  it  takes  approximately  15  hours  to  finish  the 
complete  flux-linkage  'mapping'  of  a  6/4  SRM.  If  a  two-phase  excitation  is  assumed,  the  number 
of  required  FEA  points  in  order  to  have  the  same  accuracy  will  be  3375  (15  different  rotor 
positions  times  15  different  ii  levels  times  15  different  h  levels)  requiring  225  hours  of  simulation 
(approximately  10  days,  24  hours  a  day)  and  a  large  amount  of  storage  space.  Therefore,  a 
different  approach  is  required  for  multi-phase  modeling  of  the  SRM  to  reduce  the  number  of 
needed  data  points. 


Figure  A-5  The  6/4  SRM  flux-linkage  characteristic  M  Q,  i). 


The  modeling  of  a  single-phase  excited  SRM  is  based  on  the  torque  equation  (i.e.,  (A-3) 
and  (A-4)),  and  the  voltage  equation  of  the  stator  phase,  which  is  expressed  as: 


±V(t)  =  Ri(t)  + 

dt 


(A-6) 


where  V  is  the  voltage  applied  to  a  phase  and  R  is  the  phase  resistance. 

Several  characteristics  can  be  used  to  classify  single-phase  excited  SRM  modeling 
techniques.  Here,  the  approach  solving  (A-6)  is  utilized  as  a  classification  parameter.  As  Figure 
A-6  illustrates,  three  different  categories  are  determined,  namely:  the  flux-linkage  approach,  the 
inductance  approach,  and  the  phase-current  approach. 

From  Figure  A-6,  it  is  noted  that  the  flux-linkage  approach  is  based  on  one  of  the  most  general 
solutions  of  (A-6)  which  is  the  calculation  of  a  flux-linkage  equation  M6(t),i(t))  to  be  derived 
with  respect  to  time.  This  approach  is  simple;  however,  it  is  very  sensitive  to  measurements  or 
calculations  errors  of  the  data  points  since  the  solution  is  based  on  the  derivative  of  a  function 
obtained  through  curve-fitting  techniques. 

A  different  solution  is  adopted  on  the  inductance  approach  where  equations  for  the 
derivative  of  the  flux  linkage  with  respect  to  rotor  position  (back-EMF)  and  derivative  of  the  flux 
linkage  with  respect  to  current  (inductance)  are  obtained.  This  solution  avoids  the  sensitivity 
problem  of  the  previous  solution;  however,  it  requires  solving  two  equations  instead  one. 

Finally,  in  the  phase-current  approach  the  flux  linkage  is  obtained  through  integration  requiring 
an  expression  for  the  phase  current  i(6(t),M0)  as  shown  in  Figure  A-6.  This  approach  avoids  the 
sensitivity  problem  of  the  first  approach  and  presents  a  simpler  solution  than  the  second 
approach. 
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Figure  A-6  Classification  of  modeling  techniques  for  single-phase  excited  SRM. 


In  the  next  three  sub-sections,  the  three  approaches  (flux-linkage,  inductance  and  phase- 
current)  are  explored  presenting  an  historical  evolution  of  each  approach,  a  description  of  the 
different  used  curve-fitting  techniques  and  the  data  .requirements. 


A.3.1  The  flux-linkage  approach 

A  mathematical  expression  for  the  flux  linkage  A.(9(t),i(t))  is  first  obtained  using  different 
curve-fitting  techniques.  Next,  the  derivative  with  respect  to  time  of  the  flux 

linkage  dt  is  obtained  in  order  to  solve  (A-6). 

Several  curve-fitting  techniques  or  combination  of  these  have  been  used  to  approximate 
the  flux  linkage  M9(t),i(t)),  namely:  Fourier  series,  exponential  series,  polynomial,  etc.  Several 
approaches  were  summarized  by  Fulton  and  Stephenson  in  [A.  1  ].  Here,  the  two  principal 
techniques  are  considered:  namely 

•  Polynomial  and  Fourier  series  approximations. 

•  Sum  of  exponential  approximations. 


A.3.1.J  Polynomial  and  Fourier  series  approximations 

Pickup  and  Tipping  presented  one  of  the  first  descriptions  of  the  flux  linkage  M9,i)  using 
polynomial  approximations  [A.2],  In  this  work,  both  the  rotor  position  and  phase  current  are 
approximated  by  a  sum  of  exponential  functions.  In  a  latter  method,  a  Fourier  series  is  used  to 
approximate  the  flux  linkage  X(9,i)  as  a  function  of  the  rotor  position,  and  a  polynomial 
approximation  to  describe  the  flux  linkage  X(9,i)  as  a  function  of  the  phase  current.  Analytical 
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30  di  30  J 

obtained  analytical  expression  for  M9,i).  However,  the  time-consuming  process  of  obtaining  the 
complete  description  of  the  flux  linkage  vs.  current  with  respect  to  rotor  position  is  required  to 
obtain  the  coefficients  of  these  analytical  expressions.  This  'mapping'  process  consists  on  the 
measuring/calculation  of  the  flux  linkage  (z-axes)  at  ^-degree  rotor-position  intervals  and  y-A 
winding-current  intervals  (see  Figure  A-5). 

A3. 1.2  Sum  of  exponential  approximations 

This  approach  is  presented  by  Ilic-Spomg  et.  al  in  [A.3].  Lately,  Byne  and  Dwyer 
presented  in  [A.4]  a  simplified  version  of  the  same  approach.  In  this  work,  the  sum  of  three 
exponential  functions  is  used  to  fit  the  curve  M9,i)\g=cte  using  the  flux  linkage  measured  at 
equally-spaced  rotor  positions  and  phase  currents. 
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determined  by  means  of  linear  interpolation.  However,  the  complete  description  of  the  flux 
linkage  vs.  current  with  respect  to  rotor  position  is  required  as  in  A.  3. 1.1. 


The  exponential-sum  concept  with  some  slightly  modifications  is  used  by  Torrey  et.  al  in 
different  publications.  In  [A.5],  the  exponential-sum  concept  is  incorporated  in  the  simulation  of  a 
60  kW  SRM.  However,  once  again,  the  complete  description  of  the  flux  linkage  vs.  current  with 
respect  to  rotor  position  is  required  to  obtain  the  coefficients  of  the  exponential  functions.  In 
addition,  the  SRM  is  only  simulated  at  constant  speed.  In  a  later  publication  [A.6],  relations 
between  the  coefficients  of  the  exponential  sum  and  the  physical  structure  of  the  SRM  are 
presented.  The  goal  of  that  work  is  the  reduction  of  the  number  of  measured  data  or  the  need  for 
FEA  simulations.  However,  an  iterative  analysis  is  required  to  obtain  certain  parameters. 
Moreover,  there  are  differences  between  measured  and  predicted  waveforms. 


In  summary,  the  latest  version  of  the  flux-linkage  approach  reduces  successfully  the 
number  of  required  data  points;  however,  the  sensitivity  of  this  method  to  measurements  or 
calculations  of  the  data  point  errors  is  still  present.  In  order  to  eliminate  this  sensitivity  problem, 
curve-fitting  techniques  are  used  to  approximate  the  two  partial  derivatives  of  the  flux  linkage  in 
the  inductance  approach. 


A.3.2  The  inductance  approach 


The  two  partial  derivatives  of  the  flux  linkage,  the  inductance —{e,i) 
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presented  in  [A.7]  a  method  based  in  the  same  principle  used  by  Pickup  and  Tipping  in  [A.2].  In 
this  method,  a  Fourier  series  as  a  function  of  rotor  position  and  a  polynomial  approximation  as  a 
function  of  phase  current  are  used  to  approximate  the  inductance  and  the  back-EMF,  respectively. 
The  coefficients  of  the  Fourier  series  are  obtained  from  a  relation  among  the  values  of  the 
inductance  at  three  different  rotor  positions:  aligned,  unaligned  and  midway.  The  coefficients  of 
the  polynomial  approximation  used  to  represent  the  relationship  between  phase  inductance  and 
the  back-EMF  are  derived  using  few  data  points  from  FEA  or  measurements.  Therefore,  the 
complete  description  of  the  flux  linkage  vs.  current  with  respect  to  rotor  position  is  not  required. 
Finally,  an  analytical  expression  for  torque  T(6,i)  is  determined  using  the  inductance  and  back- 
EMF  analytical  expressions. 


The  principal  contribution  of  this  approach  is  the  elimination  of  the  data-error  sensitivity 
problem  of  the  flux-linkage  approach.  In  addition,  the  amount  of  data  required  from  FEA 


simulations  or  measurements  is  reduced.  However,  this  approach  presents  a  more  difficult 
implementation  and  a  smaller  accuracy  than  the  flux-linkage  approach. 

In  order  to  combine  the  advantages  of  the  flux-linkage  and  inductance  approaches  (A-3) 
is  solved  using  a  complete  different  approach  which  is  introduced  in  the  next  section. 


A.3.3  The  phase-current  approach 

The  phase  voltage  equation,  (A-3),  is  reformulated  as  follows: 


dM  B,i ) 

dt 


=  ±V-R*i(6,i) 


(A-7) 


This  non-linear  differential  equation  is  solved  through  numerical  integration;  instead  of 
derivation  as  it  was  done  in  the  two  previous  approaches.  This  improves  the  solution  accuracy 
and  robustness  to  numerical  problems. 

Several  techniques  have  been  used  to  approximate  the  phase  current  required  in  the 
integration  of  (A-7).  Here,  the  three  principal  techniques  are  only  considered:  namely: 

•  Non-linear  differential  equation 

•  Piecewise  low-order  polynomial 

•  Gage  curve 


A.  3. 3.1  Non-linear  differential  equation 

The  idea  of  presenting  (A-6)  as  a  non-linear  differential  equation  is  initially  presented  by 
Stephenson  and  Corda  in  [A.8],  The  main  objective  of  this  approach  is  to  avoid  the  fundamental 
difficulty  of  obtaining  accurate  values  for  the  differential  of  a  function  that  is  only  specified  by 
tabulated  points.  The  solution  of  (A-7)  is  given  by: 

t 

X(d(t),i(t))=  J  [±V(t )-Ri{6(t)Mt ))]dt  (A-8) 

0 

At  each  step  of  the  numerical  integration  of  (A-8),  a  value  of  i(X,Q)  is  calculated 
performing  a  quadratic  interpolation  in  A,  and  a  linear  interpolation  in  6  from  stored  values. 


The  coenergy  is  calculated  using  numerical  integration  of  W'(6,i )  =  jA((9,/).J/|e=conjt , 

o 


and  subsequently  the  torque  T(0,i)  = 


dW'(0,i ) 


dB 


is  determined  from  the  expression  of 


coenergy  through  numerical  differentiation.  This  method  avoids  the  calculation  of  differential 


coefficients.  However,  the  complete  description  of  the  flux  linkage  vs.  current  with  respect  to 
rotor  position  is  required.  In  addition,  the  accuracy  of  the  torque  waveform  is  reduced  because  of 
the  numerical  differentiation. 

A  3. 3. 2  Piecewise  low-order  polynomial 

The  piecewise  low-order  polynomial  concept  is  presented  by  Miller  in  [A.91.  Here,  (A-7) 
is  solved  using  numerical  integration  as  in  the  previous  method;  however,  instead  of  a  set  of 
stored  values,  a  'sinusoidal'  interpolation  is  used  to  calculate  the  current  value  of  i(X,  9)  required  at 
each  step  of  the  numerical  integration  as  shown  in  Figure  A-7.  Therefore,  the  complete 
description  of  the  flux  linkage  vs.  current  with  respect  to  rotor  position  is  not  required  any  longer. 
However,  once  again,  the  torque  waveform  is  obtained  through  the  numerical  differentiation 
reducing  its  accuracy. 

In  order  to  improve  the  accuracy  of  the  torque  waveform,  an  analytical  expression  instead 
a  numerical  differentiation  is  required.  The  concept  of  the  gage  curve  introduced  in  the  next 
section  allows  obtaining  such  an  analytical  torque  expression. 

A.3.3.3  Gage  curve 

The  use  of  this  concept  in  SRM  applications  is  introduced  by  Miller  in  [A.  10].  Gage 
curves  are  used  to  get  the  instantaneous  current  values  required  on  the  calculation  of  (A-5)  from 
the  known  rotor  position  and  flux  linkage. 

As  shown  in  Figure  A-8a,  this  concept  is  similar  to  the  sinusoidal-interpolation  concept. 
However,  the  current  level  is  used  as  a  parameter  (instead  of  the  flux  linkage)  as  shown  in  Figure 
A-8b. 


Sinusoidal 

interpolation 


Unaligned 


Figure  A-7 a  Flux-linkage  X  vs.  current  with  9 as  parameter. 
A-7b  Sinusoidal  interpolation  for  phase  current  i. 


Figure  A-8a  Flivc-linkage  X  vs.  current  with  9  as  parameter.  A-8b  Gage  curve  for  different 

current  levels. 


As  in  the  sinusoidal-interpolation  concept,  mathematical  expressions  for  the  unaligned 
and  aligned  flux-linkage  curves  are  required  in  order  to  define  the  gage  curve.  However,  using  the 
gage  curve,  the  original  flux  linkage  can  be  expressed  as  the  product  of  two  functions.  Thus,  the 
integration  of  the  flux  linkage  (A-4)  yields  an  explicit  expression  for  the  coenergy  allowing 
deriving  an  expression  for  the  instantaneous  torque.  The  principal  advantage  of  this  technique  is 
that  phase  current  and  torque  are  obtained  from  an  analytical  equation  eliminating  numerical 
differentiation  and  the  need  for  a  large  number  of  FEA  simulations.  In  addition,  a  small 
computational  time  is  achieved  because  of  the  simple  formulation  of  the  gage-curve  equations. 
The  disadvantage  of  this  method  is  the  empirical  formulation  of  the  gage  curve.  Consequently, 
there  is  not  a  strong  relation  between  the  equations  and  the  SRM  physical  parameters. 

The  gage-curve  concept  is  also  used  in  [A.  11]  where  the  set  of  required  data  is  reduced, 
and  a  new  method  to  model  the  unaligned  and  aligned  flux-linkage  curves  is  presented.  Finally, 
an  ultra-fast  SRM  model  based  on  the  gage-curve  concept  is  presented  in  [A.  12]  where  the  gage 
curve  is  normalized  to  allow  a  simpler  formulation  and  a  smaller  computational  time. 

Based  on  its  numerous  advantages,  the  gage-curve  method  was  initially  selected  in  this 
Research  to  model  a  6/4  SRM.  However,  implementing  the  ultra-fast  gage-curve  modeling 
technique  [A.  12]  in  Matlab/Simulink  requires  the  solution  of  complex  mathematical  expressions 
demanding  the  use  of  different  Matlab  functions;  therefore,  the  Simulink  Matlab  Function  block 
has  to  be  repeatedly  used.  In  order  to  perform  the  different  calculations,  the  Simulink  Matlab 
Function  block  sent  the  function  inputs  data  from  the  Simulink  workspace  to  the  Matlab 
workspace  where  the  function  is  executed  as  any  Matlab  function.  Once  it  is  executed,  the  outputs 
are  saved  on  the  Matlab  workspace  and  sent  back  to  the  Simulink  workspace  to  be  the  outputs  of 
the  Matlab  Function  block.  This  exchange  of  information  between  Simulink  and  Matlab 
workspaces  has  a  negative  impact  in  the  simulation  computational  time.  Moreover,  time- 
consuming  functions  such  as  cubic-spline  interpolation,  polynomial  fitting,  or  polynomial-roots 
solution  are  normally  required;  therefore,  the  computational  time  is  further  increased.  A  new 
modeling  technique  is  introduced  in  the  next  section  in  order  to  avoid  the  above-mentioned 
disadvantages. 


This  modeling  technique  solves  (A-7)  through  numerical  integration;  therefore,  it 
classifies  under  the  current  approach  category.  However,  the  'look-up  tables'  modeling  technique 
is  classified  here  under  a  new  category  since  this  technique  will  be  used  as  comparison  through 
this  entire  report. 

The  operation  principle  of  this  technique  is  very  simple;  (A-6)  is  solved  through 
numerical  integration  to  obtain  the  flux  linkage  as  done  in  any  phase-current  approach.  However, 
the  required  values  of  i(X,6)  are  now  calculated  using  a  Simulink  look-up  table  which  performs  a 
cubic  interpolation  between  the  stored  data.  The  flux  linkage  and  the  rotor  position  are  required  as 
inputs. 


Another  Simulink  look-up  table  is  used  to  obtain  the  instantaneous  torque  T(i,  9)  in  which 
the  phase  current  and  rotor  position  are  required  as  inputs.  Figure  A-9  shows  the  Simulink 
schematic  implementation  of  the  'look-up-table'  based  model. 

The  modeling  technique  based  on  'look-up  tables'  is  able  to  overcome  almost  all  the 
disadvantages  of  previous  methods  presenting  a  very  simple  structure  and  implementation  and 
high  accuracy.  In  addition,  this  technique  provides  an  excellent  link  between  static  (FEA)  and 
dynamic  simulations.  Finally,  all  mathematical  expressions  are  implemented  using  'pure' 
Simulink  blocks  improving  further  the  computational  time. 
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Figure  A-9  Schematic  block  diagram  of  the  'look-up  table'  modeling  technique. 


The  disadvantage  of  this  modeling  technique  is  that  the  complete  description  of  (a) 
current  vs.  flux  linkage  with  respect  to  rotor  position,  and  (b)  torque  vs.  current  with  respect  to 
rotor  position  are  required.  Finally,  note  that  even  thought  the  full  'mapping'  of  the  SRM  flux 
linkage  by  means  of  FEA  is  a  time-consuming  process,  the  time  required  to  obtain  the  data  (a) 
using  FEA  is  the  same  as  the  time  required  to  obtain  the  data  (a)  and  (b)  since  both  results  can  be 
obtained  at  the  same  time.  In  addition,  it  is  assumed  that  the  SRM  physical  dimensions  are 
designed  using  FEA  results;  therefore,  the  code  required  for  the  FEA  simulation  is  available 
eliminating  one  of  the  most  time-consuming  parts  of  a  FEA  simulation. 


The  application  area  and  characteristics  of  the  SRM  model  should  be  established  in  order 
to  select  the  dynamic  modeling  technique.  FEA-based  software  packages  has  been  proving  to  be 
the  most  suitable  type  of  software  for  SRM  static  simulations.  FEA-based  software  are  used  for 
the  design  of  the  motor  itself  and  the  analysis  of  different  SRM  topologies  such  as  short-,  full- 
and  fractional-pitch  windings,  multi/single-phase  excitation,  etc.  On  the  other  hand,  circuit-based 
software  packages  such  as  SPICE  or  Saber  are  more  suitable  for  the  design  of  the  converter. 
Using  this  software,  different  converter  topologies,  snubber,  gate  and  protection  circuits  are 
studied.  Finally,  a  SRM  model  capable  of  accurate  prediction  of  the  main  SRM  waveforms  and  a 
small  computational-time  is  required  in  the  design  of  the  controller  and  control  algorithm. 

In  a  traditional  SRM  design  procedure,  a  preliminary  design  is  obtained  using 
mathematical  expressions  and  static  FEA-based  simulations.  All  the  SRM  physical  dimensions  as 
well  as  parameters  such  as  static  torque,  number  of  turns,  current  and  voltage  ratings  are  obtained 
on  this  first  step  of  the  design  procedure.  If  such  FEA-based  simulations  were  used,  the  remaining 
application  area  of  the  SRM  model  is  the  prediction  of  the  dynamic  behavior  to  develop  the 
controller  and  control  algorithms.  Therefore,  the  dynamic  modeling  technique  selected  should 
have  the  following  characteristics: 

•  Accurate  prediction  of  T(t),  X (t),  i(t),  aft),  and  6(t). 

•  Capability  of  modeling  different  SRM  topologies  (6/4,  8/6,  etc.)  as  well  as 
different  types  of  winding  pitch;  i.e.,  short,  full  and  fractional. 

•  Capability  of  modeling  different  types  of  excitation;  i.e.,  single-  and  multi¬ 
phase  excitation. 

•  Small  computational-time  requirements  to  reduce  the  developing  time  of  the 
control  algorithm. 

•  Interaction  with  FEA-based  software. 

•  Reduced  setup  time. 

•  Reduced  set  of  parameters  and  data  points. 


•  Flexibility  for  future  modifications  or  additions;  i.e.,  loss  calculations. 

The  modeling  technique  based  on  the  gage-curve  concept  is  one  of  the  most  suitable 
techniques  to  accomplish  some  of  the  above-mentioned  features.  The  particularly  simple  structure 
of  the  gage  curve  is  the  key  to  not  only  the  modeling  of  the  main  non-linearities  related  with  the 
variable  geometry  of  the  SRM  but  also  to  small  computational-time  requirements.  However,  the 
gage  curves  are  based  on  empirical  equations;  therefore,  the  design  of  the  SRM  using  this 
technique  is  not  straight  forward  since  the  results  of  the  simulation  cannot  be  easily  translated 
into  the  physical  parameters  of  the  SRM.  In  addition,  the  gage-curve  modeling  technique  has  a 
complex  setup;  consequently,  its  understanding  and  implementation  are  time  consuming. 
Moreover,  the  advantage  of  the  reduced  data  set  required  from  FEA  simulations  is  reduced  by  the 
fact  that  one  of  the  most  difficult  steps  in  a  FEA-based  simulation  is  its  setup  (code  writing  and 
debugging).  In  other  words,  the  man-hours  required  for  the  setup  of  a  few  FEA  simulation  points 
or  some  hundred  points  are  almost  the  same. 

To  overcome  the  disadvantages  presented  by  the  gage-curve  modeling  technique,  the 
'look-up-table'  modeling  technique  is  developed  and  implemented  in  this  research.  Among  the 
most  outstanding  advantages  of  this  modeling  technique  are  its  extremely  simple  implementation 
and  its  improved  computational  efficiency  compared  with  the  gage-curve  modeling  technique. 
The  main  disadvantage  is  the  larger  number  of  data  required  from  FEA  simulations;  however,  this 
data  is  available  since  the  SRM  is  assumed  to  be  designed  using  FEA.  Therefore,  the  benefits  of 
the  'look-up-table'  modeling  technique  outnumber  its  drawbacks.  The  last  step  on  the  selection  of 
the  SRM  modeling  technique  is  the  selection  of  the  software  used  to  implement  the  two-selected 
modeling  technique.  This  is  done  in  the  next  section. 


The  software  selected  to  implement  the  SRM  dynamic  model  should  allow  to  handle 
several  parts  of  the  system  that  are  very  different  in  nature.  In  general,  converters  are  modeled 
using  its  circuit  representation.  Controller  and  control  algorithms  are  expressed  in  terms  of 
mathematical  expressions  such  as  differential  or  difference  equations,  transfer  functions,  etc. 
Mechanical  loads  are  modeled  using  differential  equations.  Finally,  mathematical  tools  such  as 
the  Laplace  transform,  interpolation,  integration,  differentiation  will  be  required  for  the  modeling 
of  an  SRM  drive.  Several  simulations  package  such  Saber  in  [A.  13]  and  Pspice  in  [A.14-A.15] 
have  been  used  on  the  implementation  of  the  dynamic  simulation  technique.  However, 
Matlab/Simulink  is  selected  here  to  implement  the  two  modeling  techniques,  gage  curves  and 
'look-up'  tables,  because  Matlab/Simulink  meets  all  the  above  mentioned  conditions  and  presents 
the  following  characteristics: 

•  Powerful  tool  for  analyzing  and  developing  the  controller  and  control 
algorithms. 

•  Programming  capability  (Matlab,  C)  with  extended  libraries. 

•  Interface  with  other  programs  that  can  be  used  in  the  development  of  other 
parts  of  the  system,  such  as  Saber  in  the  development  of  the  converter. 


•  Large  Mathematical  libraries. 

•  Worldwide  known. 

Finally,  other  modeling  techniques  to  be  developed  in  the  future  can  be  easily 
implemented  in  Matlab/Simulink. 


In  this  Section,  the  basic  mathematical  description  of  a  single-phase  excited  SRM  and  the 
difficulty  of  an  accurate  modeling  of  the  SRM  were  firstly  presented  introducing  the  need  of  an 
integrated  CAD  tool  in  the  design  of  a  SRM  drive.  Next,  a  classification  of  different  modeling 
approaches  was  performed.  Table  A-l  summarizes  the  advantages  and  disadvantages  of  the 
modeling  approaches  reviewed  in  this  Section.  Two  modeling  techniques,  the  'gage-curve'  and 
'look-up-tables',  were  selected  to  be  implemented  in  Section  B.  Finally,  Matlab/Simulink  was 
selected  as  the  software  package  to  be  used  in  the  implementation  of  the  selected  simulation 
techniques. 


Table  A-l  Advantages  and  disadvantages  of  different  SRM  modeling  techniques. 


Modeling  Technique 

Advantages 

Disadvantages 

Flux  Linkage 

Polynomial  and  Fourier  series 
approximation 

Sum  of  exponential 
approximation 

Simplicity 

Large  number  of  input  data 

Limited  accuracy  because  of 
numerical  differentiation 

Large  computational-time 
requirements 

Inductance 

Reduced  number  of  input  data 

Elimination  of  numerical 
differentiation 

Complex  mathematical 
description 

Limited  accuracy  because  of 
pour  behavior  of  the  curve 
fitting 

Phase  Current 

Non-linear  differential 
equation 

Piecewise  low-order 
polynomial 

Elimination  of  numerical 
differentiation 

Large  computadonal-time 
requirements  because  of 
computing  the  integral  of 
coenergy 

Gage  curve 

Small  computational-time 
requirements 

Small  number  of  input  data 

Complex  mathematical 
description  and  understanding 

Long  setup  time 

Empirical  method 

Limited  accuracy 

Look-up  table 

Simplicity 

Accuracy 

Very  small  computational¬ 
time  requirements 

Very  small  setup  time 

Large  number  of  input  data 
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B  IMPLEMENTATION  OF  MODELING  TECHNIQUES  FOR  SINGLE-PHASE-EXCITED 
SRM 


Several  modeling  techniques  predicting  the  dynamic  behavior  of  a  single-phase-excited 
SRM  were  introduced  in  Section  A.  Each  modeling  technique  has  different  application  areas, 
accuracy,  computational-time  requirements  and  ability  to  model  non-linearities  associated  with 
the  SRM  structure.  Two  different  modeling  techniques  were  selected;  namely,  the  ’gage-curve’ 
and  look-up-table’  modeling  techniques. 

The  modeling  technique  based  on  the  gage-curve  concept  was  selected  because  of  its 
right  balance  among  data  requirements,  computational-time  requirements,  and  accuracy  on  the 
prediction  of  the  principal  waveforms  (flux  linkage,  phase  current,  torque,  etc). 

The  modeling  technique  based  on  the  look-up  tables’  concept  was  selected  because  it 
presents  a  very  simple  structure,  small  computational  time,  and  high  accuracy  on  the  prediction  of 
the  principal  waveforms.  In  addition,  this  modeling  technique  provides  an  excellent  link  between 
static  (FEA)  and  dynamic  simulation  results  as  it  will  be  shown  in  this  Section. 


The  fundamental  relationships  among  flux  linkage,  phase  current  and  rotor  position  are 
given  by  the  magnetization  curves  as  shown  in  Figure  B-la.  The  gage  curves  can  be  seemed  as  a 
graphical  transformation  of  the  classical  magnetization  curve.  The  gage  curves  are  plotted  as  flux 
linkage  vs.  rotor  position  with  current  as  parameter;  its  representation  lies  in  a  normal  plane  with 
respect  to  the  magnetization  curve  plane  [B.l].  Using  this  new  representation,  the  original  flux 
linkage  can  be  expressed  as  the  product  of  two  functions,  M.h&)  =  fl(i)  f2(6)  which  can  be 
incorporated  into  (A-4)  to  be  integrated.  Therefore,  it  is  not  necessary  to  evaluate  the  integral  of 
the  coenergy  (A-4)  at  every  time  step  of  the  simulation  process  since  an  explicit  expression  for 
the  instantaneous  torque  can  be  derived  from  the  coenergy  expression  [B.2]. 

As  shown  in  Figure  B-lb,  the  gage  curve  changes  its  size  and  shape  with  phase  current 
level  (or  saturation  level).  Thus,  the  gage  curve  should  model  accurately  these  changes  in  shape 
and  remain  as  mathematically  simple  as  possible;  this  is  the  essential  requirement  for  a  small 
computational  time. 


Figure  B-la  Flux  linkage  A  vs.  phase  current  i  with  rotor  position  das  parameter.  B-lb  Gage 

curve  for  two  different  phase  current  levels. 

B.2.1  Incorporation  of  the  Gage-Curve  Concept  Into  Dynamic  Simulations 

In  the  gage-curve  modeling  technique,  (B-3)  that  is  repeated  here  for  convenience  as  (B- 


±V  =  Ri  + 


dMdj} 

dt 


is  solved  through  numerical  integration  in  order  to  obtain 


X (0,i)  =  j  [±V -Ri( 6, A )]dt 

0 


(B-l) 


(B-2) 


where  V  is  the  applied  phase  voltage,  i  is  the  phase  current,  A  is  the  flux  linkage  and  6  is  the  rotor 
position. 

The  phase-current  value  that  is  required  at  each  time  step  on  the  numerical  solution  of  (B- 
2)  is  determined  from  the  gage  curve  using  ’past’  values  of  the  flux  linkage  and  rotor  position.  In 
addition,  mathematical  expressions  are  derived  to  calculate  the  instantaneous  torque.  Finally, 
speed  and  position  are  calculated  using  the  following  mechanical  equations: 


TJt  )—Tl(  t  )  =  Jt 


d2d  d6 

— t—  +  Bt - 

dt2  T  dt 


(B-3) 


0 


(B-4) 


dt 


where  TL  is  the  load  torque,  JT  is  the  total  inertia  of  the  system,  BT  is  the  total  friction  coefficient, 
co  is  the  instantaneous  SRM  rotor  speed,  and  do  is  an  arbitrary  initial  rotor  position. 


B.2.2  Mathematical  Formulation  of  the  Gage  Curve 

As  shown  in  Figure  B-lb,  gage  curves  change  in  size  and  shape  with  phase  current  level. 
The  gage  curves  are  normalized  in  order  to  be  independent  of  size  changes  and  simplify  the 
mathematical  definition  and  solution  of  the  equations.  This  normalization  is  made  in  such  a  way 
that  if  there  were  no  magnetic  saturation,  there  would  be  only  one  normalized  curve  for  all 
current  levels  [B.3];  meaning  that  the  gage  curve  would  be,  in  addition,  independent  of  shape 
changes.  Unfortunately,  SRM  works  under  deep  saturation,  requiring  that  the  normalized  gage 
curve  models  those  shapes  changes  accurately  in  order  to  maintain  the  accuracy  of  the  modeling 
technique. 

The  flux  linkage,  A,  is  normalized  as  follows: 

y(x)  =  ^LA  (B-5) 


where  y(x )  is  the  normalized  flux  linkage,  Aa  is  the  flux  linkage  at  the  aligned  rotor  position,  A,  is 
the  flux  linkage  at  the  unaligned  rotor  position,  and  A(x)  is  the  flux  linkage  at  any  intermediate 
rotor  position. 

The  rotor  position,  0,  is  normalized  as  follows: 


x  - 


(B-6) 


where  0a  is  the  aligned  rotor  position,  du  is  the  unaligned  rotor  position,  and  6  is  any  intermediate 
rotor  position. 

The  normalized  gage  curve  is  divided  into  five  subintervals  in  order  to  define 
representative  mathematical  expressions  for  y(x)  in  each  subinterval  as  shown  in  Figure  B-2.  The 
points  J  and  K  are  determined  by  the  pole  arc;  J  is  the  point  in  which  the  stator  and  rotor  poles 
start  to  overlap,  and  K  is  the  point  in  which  the  stator  and  rotor  poles  are  completely  overlapped 
[B-4].  Note  that  this  report  follows  the  notation  used  by  Miller  on  [B-4]. 


Region  d  extends  from  the  point  U  to  x  =  Jo  =  xJO  where  y(x)  is  defined  as: 

y(x)  =  0 

Region  e  extends  from  x  =  K)  =  xKl  to  point  A  where  y(x)  is  defined  as: 


y(x)  =  l 


Figure  B-2.  Normalized  gage  curve. 


Region  j  extends  from  x  =  JO  =  x10  to  x  =  Ji  =  xn  where  y(x)  is  defined  as: 

sin(qu ) 


/  i  mp 

y<x>=Y 


u -- 


(B-7) 


where  u  =  x-x10  and  q  =  % lj.  In  this  region,  the  first  and  second  derivatives  of  the  normalized  flux 
linkage  are  given  by: 


dy_ 

dx 


du  __  mp  r 
~dx~~T  ~ 


COS  I 


(^«)] 


mp 


qsin(  qu) 


(B-8) 

(B-9) 


Equation  (B-8)  ensures  that  there  is  no  discontinuity  in  the  rate  of  change  of  y(x)  when 


moving  to  the  adjacent  regions  d  and  c  since 


dy_ 

dx 


=  0  and 


dx 


=  m  ,  respectively. 


Region  k  extends  form  x  =  Ko  =  xKo  to  x  =  Ki  =  xKi  where  y(x)  is  defined  as: 

y(  x  )  =  l-—m0[l-sin(  rz)} 
n 


(B-10) 


where  z  =  x-Xko  and  r  =  n/2k.  In  this  region,  the  first  and  second  derivatives  of  the  normalized  flux 
linkage  are  given  by: 


dy  dy 

—  =  —  =  mn  cos(  rz) 
dx  dz  Q 

(B-ll) 

d2y  •  ,  . 

=  mQrsm(rz) 

(B-12) 

Equation  (B-ll)  ensures  that  there  is  no  discontinuity  in  the  rate  of  change  of  y(x)  when 


dy 


moving  to  the  adjacent  regions  c  and  e  since  —A  =  m0  and  —I  =  0 ,  respectively. 

dx'  1 


dy 


X=XK0 


dx 


x=xKl 


Region  c  extends  form  x  =  J  i  =  xji  to  x  =  K()  =  xKo  where  y(x)  is  defined  as: 

y(  x  )  =  av-  fiv2  +a 


(B-13) 


where  v  =  x-xj\.  In  this  region,  the  first  and  second  derivatives  of  the  normalized  flux  linkage  are 
given  by: 


-  =  —  =  a-2{3v 
dx  dv 


(B-14) 


dy 

At  P  (see Figure  B-2),  x  =  xJU  then  —  -a  =  mp  since  v  =  0. 

dx  ,_r 

X~XJ  | 


dy 

At  Q,  x  =  xK0,  then  v  =  c  and  —  -a-2 pc  =  m0 

dx  v 

X=*K0 


trip  -mQ 
2c 


The  next  step  is  to  apply  the  boundary  conditions  in  order  to  obtain  the  system  equations 
to  solve  for  the  five  unknowns  a ,  b,  h,  mP,  and  mQ  (graphically  shown  in  Figure  B-2). 


From  Figure  B-2, 


a  +  b  +  h  =  1 


(B-15) 


From  (B-7), 


(B-16) 


2 


2mnk 


From  (B-10), 

f  ^ 

<1 

II 

(B-17) 

From  (B- 13), 

j  c(mP+mQ) 

2 

(B-18) 

The  ratio  between  the  gradients  at  P  and  Q  (typically  <  1)  is  the  equation  that  defines  the 
system  completely;  that  is, 

p  =  rriQ/mp  (B-19) 

Finally,  (B-15)  can  be  rewritten  as: 


mp  = 


2 

j  +  c(l  +  M)  +  — 

n 


(B-20) 


Knowing  Jo,  Jj,  Ko,  Ki,  and  p,  the  others  values  required  for  a  complete  definition  of  the 
normalized  gage  curve  can  be  obtained  by  applying  (B-15)-(B-20). 

In  order  to  adapt  the  normalized  gage  curve  to  a  particular  SRM  topology,  FEA 
simulations  or  measurements  are  required.  As  shown  in  Figure  B-3,  the  flux  linkage  is 
measured/calculated  at  several  rotor  positions  using  two  different  current  levels:  i\  <  i,  and  12 
»is  where  is  is  the  current  at  which  the  SRM  reaches  saturation  being  at  the  aligned  rotor 
position.  These  (approximately  20)  measurements  are  used  to  plot  the  normalized  gage  curve  and 
to  calculate  graphically  the  coefficients  J0,  Jt,  Ko,  Kh  and  p.  In  addition,  the  mathematical 
formulations  of  the  unaligned  and  aligned  flux-linkage  curves  are  derived  from  the  same 
measurements.  Once  that  the  normalized  gage-curve  mathematical  structure  has  been  defined  the 
next  step  is  the  modeling  of  the  shape  variations  of  the  gage  curve  because  of  changes  on  the 
phase  current  level.  This  is  accomplished  in  the  following  section. 


B.2.3  Variation  of  the  Gage-Curve  Coefficients 

These  shape  changes  of  the  normalized  gage  curve  can  be  modeled  by  means  of  width 
changes  of  the  different  regions  (see  Figure  B-2).  At  low  phase-current  levels  ( i  <  is),  the  points  J 
and  K  can  be  assumed  fixed.  However,  the  widths  of  the  region  j  and  k  increase  with  increasing 
phase-current  level.  This  characteristic  of  the  normalized  gage  curve  can  be  taken  into  account  if 
the  points  JO  and  KO  migrate  towards  x  =  0  and  the  points  J1  and  Kl  migrate  towards  x  =  1.  The 
rate  of  variation  for  each  point  with  respect  to  the  current  level  is  assumed  constant.  This 
assumption  makes  the  model  more  empirical;  however,  it  is  considered  justifiable  because  the 
structure  of  the  gage  curve  remains  simple  and  able  to  model  the  observed  properties  of  the  SRM. 
Thus,  the  variation  of  the  breakpoints  are  expressed  as  follows: 


Figure  B-3.  Measurements  required  for  the  normalized  gage-curve  formulation. 


Jo’  —  Jo  +  Wjo  i 

For  i  >  is 

(B-21) 

Ji’  =  Ji  +Wjj  i 

For  i  >  is 

(B-22) 

Ko’  =  Ko  +  WKo  i 

For  i  >  is 

(B-23) 

K,’=  Ki  +  wK, ; 

For  i  >  is 

(B-24) 

where  WJ0,  Wji,  WKo,  and  WKi  are  the  rates  of  variation.  The  assumption  of  the  width  variation 
for  regions  j  and  k  sets  the  variation  of  the  width  for  region  c,  which  in  turn  modifies  the  value  of 
mP  as  shown  in  (B-22)  (see  Figure  B-2).  Moreover,  the  value  of  p  also  changes  from  (B-21).  The 
following  expression  considers  this  effect: 


(l-/0  (B-25) 

c 

where  %  is  the  fixed  distance  between  j  and  k,  jUo  is  the  chosen  value  of  fx  under  no  saturation 
conditions,  and  c  is  the  width  of  region  c  under  the  actual  saturation  level. 

Now  that  the  normalized  gage  curve  has  been  completely  defined,  it  will  be  used  in  the 
next  section  to  calculate  the  instantaneous  phase  current  and  torque. 


B.2.4  Instantaneous  Phase  Current  and  Torque  Calculation 

The  normalized  gage  curve  should  be  ’de-normalized’  by  means  of  mathematical 
expressions  for  the  aligned  and  unaligned  flux  linkages  in  order  to  calculate  the  instantaneous 
phase  current.  The  following  assumptions  are  made  to  define  these  mathematical  expressions  as 
illustrated  graphically  in  Figure  B-4. 


a)  The  unaligned  magnetization  curve  is  assumed  to  be  a  straight  line; 
consequently, 

Au  =  Lu  i  (B-26) 

b)  The  aligned  magnetization  curve  is  assumed  to  be  a  straight  line  for  i  <  is  and  a 
parabola  for  i  >  4;  consequently, 

Aau  =  La0  i  For  i  <  is  (B-27a) 

(h  -  As0f  =  4  a  (i  -  is0)  For  i>is  (B-27b) 

The  following  constrains  are  applied  to  ensure  a  continuous  gradient  at  i  =  is: 


Figure  B-4.  Unaligned  and  aligned  magnetization  curves. 
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(B-30) 
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where,  /Aw  =  Am  +  As,  and  -  im  -  h-  . 

andM,^;s^ 

requires  a  non-linear  calculation. 

Finally,  the  assumption  of  a  quadraUcte^  tobTnot  accurate  enough 

The  final  step  is  the  herivat.cn  C  an  equation  for 
the  flux  linkage  as  the  product  of  two  functions.  Mt.fy-fiOlfA®; 

,  .  ,  .  r  =  1J—-  from  the  integral  ot 

instantaneous  torque  can  be  derived  using  e 

coenergy,  W’=‘jA(i,0 )di .  The  normalised  gage  carve  allows  us  to  obtain  such  a  flux-linhage 
expression  (yW)  by  means  of  explicit  exp'“^  expressions 

two  regions  as  follows: 


Mx)  =  (l-y(x))  K  +  y  Aa 


If  i  <h 
From  (B-5) 

From  (B-26)  and  (B-27a) 


Mx)  =  [Lu(l-y(x))  +  y(x)Ln0]i 


(B-3 1) 


From  (A-2) 


=  Therefore, 


■  2 
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\‘M  X  )di  =J>,  (1  -  y<*»  +  1‘*  =fL- (;  •  *x))  +  y(x>  Lm  1 2 


—  (B-32) 


And 


T(6,i)=- 


dW_ 

de 


.  Consequently, 


i=const 


dW’ _dW’ dy_dx_ 

Te(6,i)~  dxdd 


(B-33) 


where, 

From  (B-32) 
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From  (B -6)  9a-6u 

Finally,  depending  upon  the  value  of  x. 
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Consequently, 


T  = 


[  La0  Lu  ]  2 


1  dy 


9„  -0U  dx 


Mx)  =  ( l-y(x))  Au  +  y(x)  Aa 


if_L>i 

From  (B-5) 

From  (B-26)  and  (B-27b) 

^x)  =  (l-y(x))Lu  i  +  y(x)[2-da^i-is0 


(B-34) 

(B-35) 


(B-36) 


(B-37) 


(B-38) 
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(B-39) 


From  (B-33) 


T.  = 


:s'  3 


J _ ty.  (B-40) 


ea-eu  dx 


Equations  (B-37)  and  (B-40)  are  the  instantaneous  torque 
smaller  and  larger  that  the  saturation  cinrem  j; "  f®  1 \  defined.  In  the  next  section,  a  ’step’ of 

- — in  °rder  to  mustrate  how  a11 

the  defined  equations  works  together. 


One  Step  of  the  Simulation  Process 


B  2*5  _ _ t. 

*  th^ion,  a  simulation  s,eP  is  solved*  "nSS 

Given  the  following  initial  conditions: 

i(t0)  -  0.  Hto)  -  0.  <W  =  4*.  **>  -  0.  VM  -  V*.  TM  -  cons, an,,  and  ,  =  0 
First  step  n_=il 

The  time  steps  of  the  SRM  dynamic  simulation  are  expressed  by: 

t,  -  to  +  n  AT  =  AT,  (B_41) 

From  (B-2)  the  new  value  of  the  flux  linkage  is  calculated  based  on  the  previous  values 
of  current  and  rotor  position  as  follows: 

(B-42) 


X(  t, )  =  }  [±  V(  t )  -  R  i(  t)]dt=  V*  A T 


The  next  step  is  to  calculate  the  ""  " 
current  is  smaller  that  the  saturation  current  or  solving  (3  38)  P 
that  the  saturation  current  is.  Consequently: 


If  i(ti)  <  is 


(B-43) 


_ Mhl _ 

1  [lu( l-y(M 1 1 )))+  y( 4h )) Ko\ 


with  y(  x(t,  ))  = 


Mh)-K 

K-K 


form  (B-5). 


If  i(ti)  >  is,  i(ti)  is  substituted  by  an  intermediate  variable  s  +<Ji  in  (B  38)  in  order  t 
express  (B-38)  as  a  quadratic  equation. 

Lu  (7 -y(x(t,)))  s2  +  2  y(x(t,))  s  +  \y(x(t,))  As0  -  MU)  +  (7-  y(x(t,)))  ho }  =  0  (B-44) 

Solving  the  quadratic  equation  in  5,  the  new  current  value  is  calculated  and  knowing  that 

i(ti)  =  s2  +  ho • 


The  value  for  the  instantaneous  torque  TM  is  calculated  applying  the present  values of 
rotor  position  and  current  to  the  instantaneous  torque  equation  derived  using  the  norma  g  g 
curve.  Consequently, 


If  i(ti)  <is 

1  dy(  x(t, ))  (B-45) 
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where>  dyjrtjjJi  will  depend  on  the  normalized  rotor  position  x(  t, )  =  -  > 
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+  *Ji(t,  )-i,  fn  -a.  -ho  f'2) 


J_ _ dy(  x(  t, ))  (B-46) 


0.  -0. 


dx 


From  (B-3),  the  speed,  aft/),  can  be  derived  as  follows: 

CO(t1)  =  j-]{Tm(t)-TL(t))dt+(Do  (B-47) 

From  (B-4)  the  new  rotor  position,  can  be  derived  as: 

0(  t1)  =  jcoIdt  +  0o  =Q)(ti  )AT  +  0O  (B-48) 

lo 


At  this  point,  all  the  unknown  values  i(t/),  co(ti)  and  6(ti)  were  calculated.  Adding 
the  known  values  V(ti),  TL(t/),  and  t2,  the  second  iteration  can  be  similarly  performed. 

Figure  B-5  shows  a  block  diagram  of  the  gage-curve-based  model  implementation  for 
only  one  phase.  For  the  multi-phase  case,  the  phases  are  considered  separately  and  the  total 
torque  is  assumed  to  be  the  sum  of  each  phase  torque  (i.e.  the  principle  of  superposition  by 
assuming  that  the  phases  are  independent). 


Figure  B-5  Block  diagram  of  the  SRM  model  based  on  the  gave-curve  concept. 


Figure  B-6  shows  the  Simulink  implementation  of  a  6/4  SRM  drive  using  the  gage-curve 
concept.  Figure  B-7  shows  the  Simulink  implementation  of  one  phase. 


Figure  B-6  Simulink  implementation  of  a  6/4SRM  drive  using  the  gage-curve  modeling 

technique. 


Figure  B-7  Simulink  implementation  of  one  phase  of  the  6/4  SRM  using  the  gage-curve 

modeling  technique. 


The  look-up  table  modeling  technique  is  based  on  simple  mathematical  expressions  and 
FEA  simulation  results.  The  FEA  results  are  incorporated  into  the  dynamic  simulation  using  look¬ 
up  tables.  The  complete  descriptions  of  the  phase  current  vs.  the  flux  linkage  with  respect  to  the 
rotor  position  as  well  as  the  torque  vs.  the  phase  current  with  respect  to  the  rotor  position  are 
required  from  FEA  or  measurements. 


The  mathematical  expressions  used  by  the  look-up  table’  modeling  technique  are  (B-2), 
(B-3),  and  (B-4).  To  obtain  the  flux  linkage,  (B-2)  is  solved  through  numerical  integration  at  each 
step  of  the  simulation  process  as  done  in  the  gage-curve  technique.  However,  the  required  values 
of  the  instantaneous  phase  current  i(t)  are  here  calculated  using  a  look-up  table’  which  performs  a 
cubic  interpolation  among  the  stored  data.  The  flux  linkage  and  the  rotor  position  are  required  as 
inputs.  Another  look-up  table  is  used  to  calculate  the  instantaneous  torque  Te(t)  in  which  the 
phase  current  and  rotor  position  are  required  as  inputs. 


Figure  B-8  shows  the  model  of  a  6/4  SRM  drive  using  the  look-up  table’  modeling 
technique.  The  multi-phase  operation  is  again  modeled  as  separated  single-phase  operation  (see 
Figure  B-9);  it  is  assumed  that  the  different  phases  do  not  have  any  influence  on  each  other.  The 
total  torque  is  then  assumed  the  sum  of  the  different  single-phase  torques. 


Figure  B-8  Simulink  implementation  of  a  6/4  SRM  drive  using  the  look-up  table1  modeling 

technique. 
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Figure  B-9  Simulink  implementation  of  one  phase  of  the  SRM  model  using  the  ’look-up 

table  ’  modeling  technique. 


Figure  B-10  shows  the  single-phase  static  torque  prediction  by  the  ’gage-curve’  and  the 
look-up  table’  models.  The  used  gage-curve  parameters  are  listed  on  Table  B-l  at  the  end  of  this 
Section.  In  order  to  have  a  closer  match  between  the  prediction  of  the  ’gage-curve’  and  the  look¬ 
up  table’  models,  the  following  gage-curve  parameters  are  modified:  Jj  =  0.4  pu  and  Lao  =  170 
mH  (instead  of  0.31  pu  and  153.1  mH,  respectively).  This  modified  gage-curve-parameter  set  is 
used  throughout  this  research.  Note  that  a  closer  match  is  achieved  by  modifying  the  gage-curve 
parameters  as  shown  in  Figure  B-10;  however,  a  difference  in  ’shape’  is  present,  especially 
between  -45°  and  -30°.  Adjusting  the  gage-curve  parameters  cannot  eliminate  this  difference  in 
shape  since  it  is  the  product  of  the  simplified  gage-curve  mathematical  definition  used  in  this 
region  [B.4]. 


Figure  B-10  Single-phase  static  torque  predicted  by  the  ’gage-curve’,  ’modified  gage-curve’ and 
look-up  table’  models  for  a  phase  current  of  13 A. 


Figure  B-ll  shows  the  dynamic  multi-phase  torque  predicted  by  the  ’modified  gage- 
curve’  and  the  look-up  table’  models.  Note  that  the  same  phase  currents  are  applied  to  both 
models;  however,  there  is  a  difference  in  the  predicted  torque.  This  difference  is  the  product  of 
the  difference  in  the  prediction  of  the  single-phase  torque  shown  in  Figure  B-10  since  both 
models  use  the  superposition  principle  to  obtain  the  multi-phase  torque. 


Rotor  position  [deg] 


Figure  B-ll  Multi-phase  dynamic  torque  predicted  by  the  ’modified  gage-curve  ’  and  look-up 

table’  models  for  a  phase  current  of  13 A. 


Figure  B-12  shows  the  flux  linkage  and  the  phase  current  predicted  by  the  ’modified 
gage-curve’  and  look-up  table’  models.  Note  that  the  flux  linkage  predicted  by  the  ’modified 
gage-curve’  model  follows  closely  the  theoretical  triangular  shape  of  the  unsaturated  flux  linkage 
(see  [B.2]  page  59)  which  is  controlled  by  the  mathematical  definition  of  the  normalized  flux 
linkage  ((B-7),  (B-10),  and  (B-13)). 


Figure  B-12  Flux  linkage  and  phase  current  predicted  by  the  ’modified  gage-curve '  and  look-up 

table  ’  models  for  a  phase  current  of  13 A. 


Finally,  note  that  the  phase  current  predicted  by  the  look-up  table’  model  decreases 
slower  than  the  one  predicted  by  the  ’gage-curve’.  It  is  because  the  aligned  inductance  ( La0 )  was 
increased  in  order  to  get  a  closer  match  of  the  torque  waveform. 


B.5  Summary  *>'  .  ;  y.  j 

Although  the  gage-curve  modeling  technique  requires  few  data  points  (typically  less  that 
20)  leading  to  a  small  computational-time  requirement,  its  implementation  and  understanding  are 
complex.  Moreover,  the  determination  of  the  used  empirical  parameters  (such  as  po ,  Wjo,  WKo)  is 
time  consuming  and  requires  previous  experience  with  the  modeling  technique.  In  other  words, 
the  parameter  adjustments  require  a  deep  knowledge  of  the  SRM  structure  because  of  the 
empirical  nature  of  the  model.  Therefore,  it  is  laborious  and  complex  for  a  person  who  did  not 
design  nor  implement  the  gage-curve  model  to  modify  it  in  order  to  accomplish  the  designer’s 
particular  needs. 


In  addition,  the  advantage  of  small  computational-time  requirements  is  reduced  on  the 
Matlab/Simulink  implementation  because  of  the  large  number  of  equations  that  has  to  be 
implemented  (more  than  25).  To  avoid  the  above-mentioned  shortcomings,  the  look-up-table’ 
based  SRM  modeling  technique  was  developed  achieving  a  simpler  to  implement  and  modify 
modeling  technique.  In  addition,  this  modeling  technique  has  smaller  computational-time 
requirements.  However,  it  requires  a  large  number  of  data  points,  typically  150  data  point,  (10 
samples  of  the  current  for  15  different  rotor  positions). 
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Table  B-l  Gage-curve  parameters  of  the  considered  6/4  SRM. 


C  MODELING  TECHNIQUES  FOR  MULTI-PHASE- 
EXCITED  SRM 


SRM  have  been  traditionally  driven  employing  single-phase  excitation;  i.e.,  only  one 
phase  is  excited  at  a  time,  except  during  the  commutation  from  one  phase  to  the  next  one.  In 
recent  years,  the  need  for  multi-phase  modeling  techniques  increased  greatly  [C.13-C-15]  since  it 
has  been  proposed  that  multi-phase  excitation  can  substantially  improves  the  performances  of 
SRM  in  terms  of  torque  ripple,  torque  density,  efficiency  and  acoustic  noise.  Moreover,  advanced 
studies  on  control  strategies  and  sensorless  operation  (torque  and  position)  increase  further  the 
need  for  a  multi-phase  model  capable  of  achieving  high  accuracy  based  on  a  reduced  number  of 
measurements  [C .  1 6-C .  24] . 

On  this  Section,  the  basic  equations  for  the  multi-phase  operation  of  the  SRM  are  first 
introduced.  There  are  two  basic  differences  between  the  single-  and  multi-phase  excited  modeling 
techniques  of  a  SRM;  a)  the  mutual-flux  interaction  among  excited  phases,  and  b)  the  total  multi¬ 
phase  torque.  When  a  single-phase  modeling  technique  is  used  to  predict  the  multi-phase 
operation  of  a  SRM,  the  mutual-flux  interaction  is  usually  neglected  and  the  total  torque  is 
assumed  to  be  the  sum  of  the  individual  torques  (i.e.,  superposition  principle)  developed  by  each 
excited  phase  [C.1-C.12].  An  example  of  this  is  the  use  of  the  'gage-curve'  or  'look-up-table' 
modeling  technique  presented  in  Section  B  as  multi-phase  modeling  techniques.  Next,  a  new 
multi-phase  modeling  technique  based  on  an  Analog  Neural  Network  (ANN)  and  ANSYS-based 
FEA  simulation  results  is  introduced.  Next,  a  6/4  ANN-based  SRM  model  is  implemented  using 
Matlab/Simulink.  Finally,  a  comparison  is  made  with  the  'gage-curve'  and  'look-up-table'  single¬ 
phase  models  implemented  in  Section  B. 


The  phase  voltage  equation  (A-3)  for  the  single-phase  excited  SRM  can  be  rewritten  for 
the  multi-phase  case  as  follows: 


dX, 

Vl=R,i,+  — 

dt 


'j  ~j 


(C-l) 


where  Vj,  Rj,  ij,  and  Xj  are  the  applied  voltage,  resistance,  current  and  flux  linkage  of  the  y-th 
phase,  respectively. 


The  flux  linkage  of  the  y'-th  phase,  A,  can  be  expressed  as  follows: 


)+£*f<hA)=y‘,.el  t+x 


jN(  i,.e,  ) 


where  Nph  is  the  number  of  phases  simultaneously  excited,  A is  the  self-flux  linkage  of  the  y'-th 
phase,  Ay  is  the  mutual-flux  linkage  between  the  y-  and  A-th  phases,  9k  is  the  relative  rotor  position 
of  the  A-th  phase  with  respect  to  the  y'-th  phase  and  A ^  represents  the  total  mutual-flux  linkage 
interacting  with  the  y'-th  phase. 

The  mutual-flux  linkage  A; n  is  a  nonlinear  function  of  the  relative  rotor  positions  and  all 
excited  currents.  This  mutual-flux  linkage  can  be  expressed  as  follows: 

=  f( )  k  *  j  (C-3) 

If  such  a  nonlinear  function  can  be  determined,  the  mutual-flux  linkage  portion  can  be 
then  subtracted  from  the  (total)  flux  linkage  in  order  to  obtain  the  self-flux  linkage  portion. 
Having  determined  the  self-flux  linkage,  single-phase-excited  SRM  modeling  techniques  can  be 
then  used  to  determine  the  phase  current.  In  other  words,  the  modeling  techniques  of  Section  A 
modified  to  account  for  the  mutual-flux  linkages  can  be  used  as  multi-phase-excited  SRM 
modeling  techniques.  Knowing  the  phase  currents,  the  total  torque  can  be  then  determined  from  a 
full  mapping  of  the  multi-phase  torque  based  on  measurements  or  FEA  results. 

Here,  the  look-up-table'  modeling  technique  is  selected  as  the  modified  single-phase- 
excited  SRM  modeling  technique  because  the  method  used  to  determine  (C-3)  is  based  on  the 
coenergy  and  field  energy  data  points  calculated  from  FEA.  Therefore,  it  is  assumed  that  the  flux- 
linkage  characteristic  for  single-phase  excitation  required  in  the  'look-up-table'  modeling 
technique  can  be  easily  obtained. 

Figure  C-l  illustrates  the  block  diagram  of  the  SRM  multi-phase-excited  model  based  on 
the  modified  'look-up-table'  in  order  to  account  for  the  phase  mutual  interaction. 

This  multi-phase  model  is  derived  from  the  single-phase  model  as  follows: 

a)  To  calculate  the  (total)  flux  linkage  of  the  y'-th  phase,  (C-l)  is  rewritten  as  follows: 


=  Vj-Rjij 


(C-4) 


b)  (C-4)  is  solved  through  numerical  integration  (as  in  any  phase-current  technique). 

c)  The  mutual-flux  linkage  is  subtracted  from  the  total-flux  linkage  to  obtain  the  self-flux 
linkage  (see  (C-2)). 


Figure  C-l  Block  diagram  for  the  SRM  multi-phase  model  based  on  the  single-phase  'look-up 

table'  modeling  technique. 


d)  The  phase  current  ij  is  calculated  using  cubic-spline  interpolation  from  values  stored  on  a 
look-up  table  that  has  the  self-flux  linkage  A#  and  the  rotor  position  6j  as  inputs  in  the 
same  fashion  as  used  on  the  'look-up-table'  modeling  technique. 

e)  The  total  torque  is  calculated  from  measurements  or  FEA  results.  It  is  assumed  that  the 
full  mapping  of  the  total  torque  is  available. 

In  summary,  the  two  basic  differences  between  the  single-  and  multi-phase  modeling 
techniques  have  been  considered  separately  by  presenting  a  new  approach  to  account  for  the 
mutual  fluxes  interactions  as  well  as  a  simplified  approach  to  obtain  the  multi-phase  torque. 
Hence,  the  problem  of  multi-phase  modeling  of  a  SRM  has  been  reduced  to  the 
derivations/calculations  of  (C-3)  and  the  multi-phase  torque.  These  issues  are  further  addressed  in 
the  next  two  sections  that  present  a  method  based  on  the  energy  and  coenergy  of  the  SRM  is 
presented. 


C.2.1  Mutual-Flux  Linkage  And  Torque  Determination 

A  multi-excited  magnetic  system  (in  this  case,  the  SRM  under  multi-phase  excitation)  is 
governed  by  the  following  relationship  among  the  flux  linkage  A„  the  phase  current  the  total 
system  coenergy  Wc,  and  the  total  system  field  energy  Wf. 


In  order  to  describe  this  method,  a  SRM  with  two  phases  simultaneously  excited  ( Nph  = 
2)  is  selected;  therefore,  (C-5)  can  be  expressed  as  [C.251: 

X,  i,  +  X2  i2  =  WC+Wf  (C-6) 

The  total-phase  fluxes  X\  and  A_>  are  divided  into  self-flux  linkages  {Xu  and  X22)  and 
mutual-flux  linkages  {Xt2  and  X2i)  as  follows: 


X  i  =  Xu  +  X/2  (C-7) 

X2  =  X22  +  A2 1  (C-8) 

Based  on  these  equations,  a  mutual  interaction  function  M  can  be  defined  as  follows: 

X,2  =  M(i,t  i2,  Q )  i2  and  X2,  =  M(U,  i2,  d )  i,  (C-9) 

Then,  (C-6)  can  be  rewritten  as: 


Xu  ii  +  X22  i2  +  2  M(ii,  i2, 6 )  it  i2=  Wc  +  Wf 


(C-10) 


The  mutual  interaction  function  M(ih  i2, 6 )  can  be  calculated  as  follows: 


M(i,,i2,0)  = 


+  W j  Xni,  X22i2 

2iih 


(C-ll) 


Note  that,  as  was  introduced  on  Section  A,  the  self-flux  linkage,  the  coenergy  and  the 
field  energy  are  function  of  the  phase  currents  and  rotor  position.  Therefore,  the  mutual 
interaction  function  M  is  a  function  of  phase  current  ii,  phase  current  i2,  and  rotor  position  6. 

Different  techniques  can  be  implemented  in  order  to  derive  the  mutual  interaction 
function  M(ii,i2l6).  One  of  the  most  simple  and  straightforward  techniques  is  the  full  mapping  of 
the  mutual  interaction  function  using  FEA  in  order  to  implement  a  tri-dimensional  look-up  table 
(for  two-phase  excitation).  However,  it  cannot  be  easily  implemented  because  of  the  large  amount 
of  required  FEA-simulations/measurements  (see  Section  A.3).  Alternatively,  a  curve-fitting 
technique  or  ANN  could  be  implemented  to  reduce  the  amount  of  information  required  from 
FEA-simulations/measurements.  Two  different  techniques  were  evaluated  in  this  report;  namely, 
Fourier  series  and  ANN. 


Fourier  series  are  a  well-known  curve-fitting  technique;  however,  it  requires  a  large 
number  of  FEA-simulation  results  (600  points  or  more)  for  this  particular  implementation;  and  its 
accuracy  for  certain  currents  levels  and  rotor  positions  is  not  satisfactory.  In  addition,  Fourier 
series  provides  a  solution  for  only  the  mutual  interaction  function  and  not  for  the  total  torque. 

ANN  is  a  new  technique  used  for  this  purpose  that  requires  particular  knowledge  of  the 
SRM  system.  Traditional  ANN  structures  can  be  trained  using  a  small  number  of  FEA-simulation 
results  (less  than  400  points)  while  achieving  satisfactory  accuracy.  However,  the  accuracy  of  the 
ANN  depends  on  several  variables  of  the  ANN  structure  (number  of  layers,  number  of  neurons, 
number  of  data  points,  type  of  transfer  function,  training  algorithm,  etc.).  Therefore,  several 


optimization  techniques  and  generic  algorithms  can  be  implemented  to  improve  the  ANN 
accuracy,  reduce  the  number  of  required  data  points,  reduce  the  time  required  for  training, 
simplify  the  ANN  structure,  etc. 

Up  to  this  point,  the  presented  multi-phase  modeling  technique  is  completely  general;  any 
SRM  topology  (i.e.,  6/4,  8/6,  10/8,  12/8,  etc.)  or  winding  configuration  (short,  full  or  fractional 
pitch)  can  be  implemented.  However,  the  SRM  topology  and  winding  configuration  modifies  the 
ANN  structures  and  the  data  points  required  for  training.  Therefore,  there  is  not  an  ANN 
structures  providing  the  best  performance  for  all  SRM  configurations;  therefore,  each 
configuration  should  be  studied  and  optimized  individually.  In  the  next  section,  the  effects  of  the 
SRM  topology  on  the  selection  of  the  ANN  structures  and  its  training  data  points  is  presented  by 
means  of  an  example. 


In  this  section,  a  multi-phase  model  of  a  6/4  SRM  is  analyzed  and  implemented  using  the 
previously  presented  multi-phase  modeling  technique.  First,  the  particular  6/4  SRM  topology  is 
analyzed  in  order  to  determine  the  structure  and  the  amount  of  data  required  to  train  the  used 
ANNs. 


The  6/4  SRM  is  essentially  designed  to  operate  under  single-phase  excitation;  its  stroke 
angle  6Sr  is  30°  and  its  positive  torque  capability  angle  is  45°  (see  Section  A.2).  Therefore, 
there  is  a  15°  region  in  which  two  phases  can  be  simultaneously  excited  producing  positive  torque 
(overlap).  As  shown  in  Figure  C-2a,  phase  1  is  the  working  phase  when  at  pair  of  rotor  poles  are 
at  -45°  (the  unaligned  position)  with  respect  to  phase  1.  At  this  point,  phase  3  (the  leading  phase) 
is  conducting;  therefore,  the  possible  interaction  is  between  the  working  phase  and  the  leading 
phase  as  the  working  phase  approaches  0°  (the  aligned  position)  as  shown  in  Figure  C-2d.  At  - 
30°,  the  leading  phase  starts  to  produce  negative  torque;  therefore,  it  is  usually  turned-off;  see 
Figure  C-2b.  At  -15°,  phase  3  is  completely  turned-off  (see  Figure  C-2c)  and  the  other  pair  of 
rotor  poles  is  at  -45°  with  respect  to  phase  2;  therefore,  phase  2  becomes  the  working  phase  and 
phase  1  the  leading  phase. 

As  shown  in  Figure  C-l  two  ANNs,  one  required  to  predict  the  mutual-flux  interaction 
(A-  =  f( , 9 6 N  ))  and  the  other  one  required  to  predict  the  torque 

(Tlola,  =  f( i, ,9,  ))>  and  a  look-up  table  are  required  to  compose  each  SRM  phase. 
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Figure  C-2a  6/4  SRM  at  -45°. 
C-2c  6/4  SRM  at -15°. 


C-2b  6/4  SRM  at  -3(f. 
C-2d  6/4  SRM  at  CT. 


In  order  to  establish  the  ANN  structures,  a  set  of  training  data  (in  this  case  calculated  by 
ANSYS-based  FEA)  is  required.  In  order  to  obtain  such  a  training  data  set,  the  current  range  is 
divide  into  5  equally-spaced  intervals  (3A,  6A,  9A,  12A  and  15A  =  Imax)  and  the  rotor-position 
range  is  divide  into  1 1  equally-spaced  rotor  positions.  Therefore,  the  number  of  used  data  points 
is  5*5*11  =  275.  Note  that  if  the  same  'density'  of  data  as  in  the  'look-up  table'  model  were 
required,  the  amount  of  required  multi-phase  data  will  be  15*15*16  =  3600  (15  current  levels  per 
phase  and  16  rotor  positions).  Finally,  the  single-phase  characteristic  used  to  implement  the  'look- 


up  table'  model  is  incorporated  to  the  training  data  set  since  it  is  already  available.  Therefore,  240 
more  points  are  added  for  a  total  of  5 15. 

In  order  to  use  just  one  software  package,  the  handling  and  training  of  the  ANN  is 
preformed  using  Matlab.  Only  two  commands  in  the  Matlab  ANN  toolbox®  are  used:  newff  and 
train.  The  newff  command  creates  a  feed-forward  ANN  while  the  train  command  trains  the  ANN 
using  the  back-propagation  algorithm.  These  commands  are  very  simple  to  use;  however,  they 
offer  a  limited  access  and  control  of  the  training  procedure. 

Several  different  combinations  of  ANN  structure,  neuron  transfer  functions  and  training 
algorithms  were  studied  for  each  type  of  ANN  (mutual  flux  and  torque).  The  selection  of  the  most 
suitable  ANN  structures  is  based  on  the  means  square  error  given  by  the  train  command.  In 
addition,  comparisons  are  made  between  a  set  of  testing  data  composed  by  35  data  points 
(obtained  in  the  same  way  as  the  training  set,  but  not  used  to  train  the  ANN)  and  the  predicted 
results  of  the  ANN.  Using  this  criteria  an  ANN  composed  of  3  neurons  in  the  input  layer,  two 
hidden  layers  with  13  and  11  neurons,  and  a  single  neuron  in  the  output  was  selected  for  the 
prediction  of  the  mutual  function  M.  For  the  prediction  of  the  torque,  the  selected  ANN  was 
composed  of  3  neurons  in  the  input  layer,  a  hidden  layer  with  13  neurons  and  a  single  neuron  in 
the  output. 

Figure  C-3  shows  the  Simulink  implementation  of  the  6/4  SRM  drive  using  the  ANN- 
based  multi-phase  technique.  Figure  C-4  shows  the  'inside'  of  the  6/4  SRM  ANN  model  of  Figure 
C-3.  In  this  Figure,  Fll,  F22  and  F33  are  the  self  fluxes  of  phases  1,  2  and  3,  respectively,  and 
FL1,  FL2  and  FL3  are  the  mutual  fluxes  calculated  by  the  ANN  allocated  inside  of  the  mutual- 
flux  calculator  block. 

The  ANN  that  predicts  the  torque  is  allocated  inside  of  the  ANN  torque  block. 


Figure  C-3  Simulink  implementation  of  the  6/4  SRM  model  using  ANNs. 


ANN  Torque 


Figure  C-4  'Inside '  of  the  SRM  block. 


Figure  C-5  shows  the  total,  mutual  and  self  fluxes  of  the  implemented  6/4  SRM  operating 
at  a  low  speed  (30  rad/sec).  The  power  converter  is  set  to  turn  the  phase  on  when  the 
'participating'  rotor-poles  pair  reaches  -45°  and  turn  the  phase  off  when  the  same  rotor-poles  pair 
reaches  -3°. 

The  single-phase  modeling  techniques,  such  as  the  'look-up  table',  neglects  the  effects  of  mutual 
interactions  among  conducting  phases  predicting  only  the  total  flux  shown  in  Figure  C-5.  This 
flux,  that  has  the  contribution  of  the  leading  phase,  is  used  to  calculate  the  working-phase  current 
using  the  'look-up  table'  block.  The  multi-phase  ANN-based  technique  substracts  the  portion 
contributed  by  the  leading  phase  (calculated  by  the  ANN)  from  the  total  flux  (see  Figure  C-4)  in 
order  to  calculate  the  flux  produced  by  the  working  phase.  Using  this  'modified'  flux,  which  is  the 
self  flux  in  Figure  C-5,  the  working-phase  current  is  calculated  using  the  'look-up  table'  block. 
Note  that  the  total  and  self  fluxes  match  when  the  leading-phase  current  is  zero. 


Figure  C-5  Total,  mutual  and  self  fluxes  of  a  phase  of  the  6/4  SRM  ANN -based  model  for  low 

speed  operation  (30  rad/sec.). 
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Figure  C-6  Total,  mutual  and  self  fluxes  of  a  phase  for  the  6/4  SRM  ANN -based  model  at  high¬ 
speed  operation  (120  rad/sec). 


Figure  C-6  shows  the  total,  mutual  and  self  fluxes  of  the  implemented  6/4  SRM  operating 
at  a  high  speed  (120  rad/sec).  The  power  converter  is  set  to  turn  the  phase  on  when  'participating' 
rotor-poles  pair  reaches  -45°  and  turn  the  phase  off  when  the  same  rotor-pole  pair  reaches  -13°. 
Note  that  both  the  total  and  self  fluxes  match  as  in  Figure  C-5  when  the  leading-phase  current  is 
zero.  It  is  because  at  high-speed  operation  (single-pulse  mode)  the  phase  current  levels  are  lower; 
therefore,  the  interaction  between  the  leading  and  working  phases  is  reduced  and  does  not  present 
an  important  effect. 
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Figure  C-7  ANN  prediction  of  the  single-phase  static  torque  and  FEA  training  data. 


Figure  C-7  shows  the  single-phase  static  torque  waveform  predicted  by  the  ANN  for 
different  current  levels  in  order  to  illustrate  the  capability  of  the  ANN-based  multi-phase  model. 
In  addition,  it  shows  the  single-phase  portion  of  the  FEA-obtained  training  data  points  used  to 
train  the  ANN.  The  single-phase  torque  is  obtained  from  the  multi-phase  trained  ANN  by  setting 
the  leading  current  to  zero.  Note  that  the  differences  among  the  ANN-predicted  waveform  and  the 
training  points  are  greater  at  low  current  levels. 


Working-phase  rotor  position  [deg] 


Figure  C-8  ANN  prediction  of  the  multi-phase  static  torque  and  FEA  training  data. 


Figure  C-8  shows  the  multi-phase  static  torque  predicted  by  the  ANN-based  model,  by 
the  'look-up  table'  model  and  the  13A  portion  of  the  ANN  training  data  set  (obtained  by  FEA). 
The  leading  and  working  phases  carry  13A  for  the  full  45°.  However,  the  difference  between  the 
'look-up  table'  model  and  the  ANN-based  model  is  smaller  at  -45°  than  at  around  -30°.  It  is 
because  the  flux  produced  by  the  working  phase  at  -45°  is  smaller  than  at  -30°  since  the  reluctance 
is  higher  at  -45°  than  at  -30°.  The  main  reason  for  this  difference  is  that  the  changes  on  the 
saturation  level  of  the  SRM  (because  of  two  phases  excited  simultaneously)  are  neglected  on  the 
single-phase  modeling  techniques  (such  as  the  'look-up  table'  modeling  technique). 
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Figure  C-9  ANN  and  'look-up  table'  predictions  of  the  multi-phase  dynamic  torque. 


Finally,  note  that  the  torque  is  equal  to  zero  at  -15°  in  both  techniques  because  at  this 
rotor  position  (see  Figure  C- 2c)  the  amount  of  attraction  produced  by  the  working  phase  (positive 
torque)  is  equal  to  the  amount  of  repulsion  produced  by  the  leading  phase  (negative  torque)  since 
both  phases  are  carrying  the  same  current. 

Figure  C-9  shows  the  multi-phase  dynamic  torque  predicted  by  the  ANN-based  model 
and  the  'look-up  table'  model.  The  dynamic  torque  predicted  by  the  ANN-based  model  presents  a 
reduction  at  certain  rotor  positions  when  compared  with  the  dynamic  torque  predicted  by  the 
'look-up  table'  model  (using  the  superposition  principle).  This  reduction  in  torque  is  produced  by 
the  increment  on  the  SRM  saturation  level  when  two  phases  are 


Figure  C-10  ANN  and  'look-up  table'  predictions  of  the  multi-phase  dynamic  torque. 


excited  simultaneously.  Note  that  the  torque  predicted  by  single-  and  multi-phase  techniques 
matches  when  the  SRM  operates  under  single-phase  excitation.  It  will  be  illustrated  later  that  this 
torque  reduction  is  also  present  in  the  measurements. 

Figure  C-10  shows  the  multi -phase  dynamic  torque  predicted  by  the  ANN-based  model 
and  the  'look-up  table'  model.  This  figure  is  similar  to  Figure  4-9.  However,  the  power  converter 
is  now  set  to  turn  the  phase  on  when  the  'participating'  rotor-pole  pair  reaches  -33°  and  turn  the 
phase  off  when  the  same  rotor-pole  pair  reaches  -13°  in  order  to  reduce  the  overlapping  between 
leading  and  working  phases.  Consequently,  the  torque  predicted  by  the  'look-up  table'  model 
(single-phase  modeling  technique)  and  the  torque  predicted  using  the  ANN-based  model  (multi¬ 
phase  modeling  technique)  match  closely  since  the  SRM  is  working  essentially  under  single¬ 
phase  excitation. 

Figure  C-ll  shows  the  multi -phase  dynamic  torque  predicted  by  the  ANN-based  model 
and  the  'look-up  table'  model  under  high-speed  operation.  Both  torque  predictions  match  closely 
for  the  same  reason  that  the  total  and  self  fluxes  match  in  Figure  C-6.  At  high-speed  operation, 
the  current  level  reached  is  lower;  therefore,  the  changes  in  the  saturation  level  of  the  SRM  are 
reduced. 
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Figure  C-ll  ANN  and  'look-up  table'  predictions  of  the  multi-phase  dynamic  torque  for  high¬ 
speed  operation. 
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Figure  C-12  Singe-phase  static  torque  comparison  between  FEA  and  measurements. 


Figure  C-12  shows  the  comparison  between  the  static  torque  predicted  by  FEA  and  the 
measured  static  torque.  It  has  been  pointed  out  that  a  2D  FEA  program  such  as  ANSYS  will 
predict  between  15%-20%  more  torque  than  the  measured  torque.  Professor  Miller  mentions  in 
his  book  Switched  Reluctance  Motors  and  their  Control  [C-l]  that:  "Two-dimensional  solvers 
cannot  be  expected  to  give  accurate  magnetization  curves  for  switched  reluctance  motors, 
especially  in  partially  aligned  position,  because  of  the  significant  influence  of  end  effects."  In 
addition,  Professor  Miller  presents  a  numerical  example  where  the  aligned  and  unaligned 
inductances  computed  using  2D  and  3D  FEA  programs  are  compared  with  measured  aligned  and 
unaligned  inductances.  The  results  of  this  computation/measurement  comparison  are  presented  in 
Table  C-l.  Note  that  the  difference  between  the  aligned  and  unaligned  inductances  using  a  2D 
FEA  is  15.3%  larger  than  difference  between  the  measured  aligned  and  unaligned  inductances.  In 
addition,  the  difference  between  the  aligned  and  unaligned  inductances  computed  using  3D  FEA 
is  only  2.8%  larger  than  the  difference  between  the  measured  aligned  and  unaligned  inductances. 
The  difference  between  the  aligned  and  unaligned  inductances  is  proportional  to  the  average 
torque  shown  in  Figure  C-12.  This  large  difference  between  2D  FEA  results  and  measurements 
strongly  suggest  that  an  accurate  multi-phase  model  of  the  SRM  should  use  3D  FEA  results. 


However,  obtaining  the  required  data  points  (around  500)  using  3D  FEA  is  impractical  since  it  is 
more  complex  and  time  consuming  than  2D  FEA.  Therefore,  measurements  are  required  as  input 
data. 


Table  C-l  Aligned  and  unaligned  inductances. 


Measured 

2D  computation 

3D  computation 

Aligned 

66  mH 

66.5  mH 

68  mH 

Unaligned 

19.8  mH 

13.2  mH 

20.5  mH 

Difference 

46.2  mH 

53.3  mH 

47.5  mH 

In  order  to  take  the  measurements,  it  is  required  to  setup  a  test/measurement  bench 
having,  in  addition  to  build  the  SRM,  a  torque  meter,  a  mechanism  that  allows  blocking  the  rotor 
at  different  positions,  two  power  supplies  with  current  limits  and  two  amperemeters.  Once  that 
the  test  bench  had  been  mounted  the  measurements  of  the  single-  and  multi-phase  flux-linkage 
characteristics  as  well  as  the  torque  characteristics  are  straightforward.  Figure  C-l 3  shows  a  s  the 
general  setup  of  the  test/measurements  bench.  Only  torque  measurements  are  presented  here;  at 
the  writing  of  this  report,  this  research  group  is  testing  a  developed  method  that  allows  measuring 
the  flux  linkage  under  multi-phase  excitation. 


Figure  C-l 3  General  setup  of  the  6/4  SRM  test/measuremnt  bench. 


Baldor  Motors  and  Drives  (Forth  Smith,  Arkansas)  built  the  6/4  SRM  used  in  the  test 


bench. 


The  equipment  used  in  the  test  bench  were  the  following: 

•  6/4  SRM. 

•  6434B  Hewlett  Packard  DC  power  supply  (MOV,  0-25A. 

•  160T60  TCR  DC  power  supply  0-120V,  0-50A. 

•  1 104-2k  Lebow  torque  sensor,  2000  lb.-in,  9000  rpm. 

•  3170  Daytronic  strain  gage  conditioner. 

•  Two  Fluke  73  multimeters. 

•  Rotor-blocking  device  capable  of  blocking  the  rotor  every  3°. 

The  torque  measurements  were  taken  in  order  to  generate  a  training  data  set  similar  to  the 
FEA-based  training  data  set.  Therefore,  the  same  phase-current  levels  were  selected  (0A,  3A,  6A, 
9 A,  12 A  and  15  A).  However,  the  measured  rotor  positions  were  every  6°  instead  of  4.5°  (0°,  3°, 
6°,  12°,  18°,  24°,  30°,  36°,  42°  and  45°);  the  total  number  of  measurements  was  6*6*10  =  360. 
Note  that  the  rotor-blocking  device  has  the  capability  of  only  blocking  the  rotor  every  3°  or 
multiple  of  3°.  Finally,  the  single-phase  torque  characteristic  was  taken  every  3°  in  order  to  have  a 
better  resolution  for  plotting  purposes;  however,  it  was  not  incorporated  into  the  measured 
training  data  set. 
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Figure  C-14  ANN  prediction  of  the  single-phase  static  torque  and  the  measured  training  data. 


Figure  C-14  shows  the  single-phase  static  torque  waveforms  for  different  current  levels 
predicted  by  the  ANN  trained  with  the  measured  training  data  set.  In  addition,  it  shows  the  single¬ 
phase  portion  of  the  measurements  used  to  train  such  an  ANN.  The  single-phase  torque  is 
obtained  from  the  multi-phase  trained  ANN  setting  the  leading-phase  current  to  zero.  Note  that, 
as  in  Figure  C-7,  the  differences  among  the  ANN  prediction  and  the  training  points  are  greater  at 
low  current  levels. 


Figure  C-15  ANN  prediction  of  the  multi-phase  static  torque  and  the  measured  training  data. 


Figure  C-15  shows  the  multi-phase  static  torque  predicted  by  the  ANN-based  model  and 
by  the  'look-up  table'  model  (by  applying  superposition  principle).  Note  that,  as  in  Figure  C-8,  the 
difference  between  the  'look-up  table'  model  and  the  ANN-based  model  is  smaller  around  45° 
(unaligned  position)  than  at  -30°.  In  addition,  the  torque  predicted  by  the  two  techniques  is  zero  at 
-15°.  Here,  the  same  explanation  presented  in  Figure  C-8  applies. 

Figure  C-16  shows  the  multi-phase  dynamic  torque  predicted  by  the  ANN-based  model 
and  the  'look-up  table’  model.  Note  that,  as  in  Figure  C-9,  the  dynamic  torque  predicted  by  the 
ANN-based  model  presents  a  reduction  at  certain  rotor  positions  compared  with  the  dynamic 
torque  predicted  by  the  'look-up  table'  model. 
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Figure  C-16  ANN  and  'look-up  table'  predictions  of  the  multi-phase  dynamic  torque. 


In  addition,  note  that  there  is  a  change  in  shape  of  the  torque  predicted  by  the  'look-up 
table'  model  on  Figure  C-15  compared  with  the  torque  predicted  by  the  'look-up  table'  model  on 
Figure  C-9.  The  reduction  of  torque  predicted  by  the  ANN-based  model  is  produced,  as  explained 
on  Figure  C-8,  by  an  increase  in  the  saturation  level  of  the  SRM  when  two  phases  conduct 
simultaneously.  The  change  in  shape  is  produced  by  the  difference  between  the  single-phase 
static  torque  predicted  by  FEA  and  the  single-phase  static  torque  measured  (see  Figure  C-12). 
The  single-phase  static  torque  measured  is  not  only  smaller  in  magnitude  but  is  also  'shifted' 
about  2°.  Since  the  torque  produced  by  the  leading  and  working  phases  is  added  in  the  'look-up 
table'  model,  this  'shift'  produces  a  peak  in  the  'look-up  table'  predicted  torque  of  Figure  C-15. 


This  Section  firstly  introduced  the  difference  between  the  SRM  single-phase  model  used 
to  represent  a  multi-phase  operation  using  the  superposition  principle  and  a  truly  SRM  multi¬ 
phase  model.  Next,  it  presented  the  basic  equations  describing  the  operation  of  the  SRM  under 
multi-phase  excitation.  Next,  it  developed  a  new  method  based  on  FEA  that  accounts  for  the 
mutual  interaction  among  conducting  phases.  This  method  accounts  for  the  changes  on  the  fluxes 
of  the  SRM  as  well  as  the  torque  when  exciting  a  second  phase.  Next,  the  developed  method  was 
implemented  using  two  ANN.  Simulation  results  (using  a  single-phase  model  -the  'look-up  table' 


model-  and  a  multi-phase  model  -the  ANN-based  model-)  were  compared  with  measurements. 
The  need  for  accurate  training  data  was  shown  based  on  those  comparisons.  Finally,  torque 
measurements  were  taken  in  order  to  generate  an  accurate  training  data  set.  Using  this  measured 
data  set,  a  new  set  of  plots  was  made  comparing  the  single-phase  model  and  the  multi-phase 
model.  The  results  of  those  comparisons  were  similar  to  the  previously  present  results;  therefore, 
the  capabilities  of  the  new  ANN-based  multi-phase  modeling  technique  presented  in  this  report  to 
model  the  behavior  of  the  SRM  under  multi-phase  operation  were  validated. 
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